Experimental investigation of classical and quantum correlations under decoherence 
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It is well known that many operations in quantum information processing depend largely on a 
special kind of quantum correlation, that is, entanglement. However, there are also quantum tasks 
that display the quantum advantage without entanglement. Distinguishing classical and quantum 
correlations in quantum systems is therefore of both fundamental and practical importance. In 
consideration of the unavoidable interaction between correlated systems and the environment, un- 
derstanding the dynamics of correlations would stimulate great interest. In this study, we investigate 
the dynamics of different kinds of bipartite correlations in an all-optical experimental setup. The 
sudden change in behaviour in the decay rates of correlations and their immunity against certain 
decoherences are shown. Moreover, quantum correlation is observed to be larger than classical 
correlation, which disproves the early conjecture that classical correlation is always greater than 
quantum correlation. Our observations may be important for quantum information processing. 

PACS numbers: 03.67. Yz, 42.50.Xa, 03.65.Ud 



I. INTRODUCTION 

Correlations, including classical and quantum parts, 
are crucial in science and technology. Although many 
operations in quantum information processing largely 
depend on a special kind of quantum correlation, that 
is, entanglement [lj, there are quantum tasks that dis- 
play the quantum advantage without entanglement 0- 
|4|, and some of these have been verified experimentally 
Q. Distinguishing classical and quantum correlations is 
therefore of both fundamental and practical importance. 
Many theoretical studies have been conducted in this di- 
rection Among them, quantifying the quantum- 
ness of correlations with quantum discord [7[ has received 
great attention (IHEIH^. 

In the field of quantum information, for a bipartite sys- 
tem Pab, it is widely accepted that quantum mutual in- 
formation measures its total correlations 0, H3| defined 
as 1{pab) = S(pa) + S(p B ) - S(pab) HI, where p A 
and pb are the reduced-density matrices of Pab- S(p) = 
— tr(plog 2 p) in the von Neumann entropy. Depending 
on the maximal information gained for pab with mea- 
surement on one of the subsystems, classical correlation 
(C) Q is defined as C(pab) = max[S(pA) - X^i^G^)] 

B]B j j 

, where B\B; is a positive-operator- valued measure per- 

formed on the subsystem B and qj = tr ab{Bj pab Bj). 

fP A = tr B(BjPABBj)/qj is the postmeasurement state 
of A after obtaining the outcome j on the particle B. 
Quantum correlation is therefore given by Q(pab) = 
I(pab) — C(pab)- In such a case, Q is just identical to 
quantum discord with a definition based on the distinc- 
tion between classical information theory and quantum 
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information theory [3], which can be further distributed 
into entanglement and the additional quantum correla- 
tion (non-entanglement quantum correlation) (29j . 



Because of the unavoidable interaction between a 
quantum system and its environment, understanding the 
dynamics of different kinds of correlations has stimu- 
lated great interest. Some studies have focused on the 
comparison between the dynamics of quantum discord 
and entanglement under both Markovian and Non- 
Mar kovian [2lL l22j environments. Their behaviours have 
been shown to be very different, and the quantum discord 
is more robust against decoherence than entanglement in 
the Markovian evolution Pt . Recently, several pecu- 
liar properties in the dynamics of classical and quantum 
correlations have also been shown with the present of 
Markovian noise [23|, 24] , in which the decay rates of cor- 
relations may exhibit sudden changes in behaviour (23j 
and the bipartite quantum correlation may completely 
disappear without being transferred to the environment 
241. 



In this study, we experimentally investigate the dy- 
namics of classical and quantum correlations between 
biqubit systems in a one-sided phase-damping channel. 
The sudden changes in behaviours in the decay rates of 
classical and quantum correlations are shown in an all- 
optical experimental setup, in which the dephasing envi- 
ronment is simulated by birefringent quartz plates. Clas- 
sical and quantum correlations are observed to remain 
unaffected under certain decoherence areas. Moreover, 
quantum correlation is shown to be larger than classical 
correlation during the dynamics of a special input state, 
which contradicts the early conjecture that classical cor- 
relation is always greater than quantum correlation [30| . 
Our observations may have important roles in quantum 
information processing. 
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II. RESULTS 

Theoretical schemes. Consider a two-level quantum 
system S with lower and upper states |0)s and \ X)s under 
the action of a phase-damping environment E (\Q)e is 
the initial state). The interaction quantum map can be 
written as [HI 



|0>s®|0)j3 -► ^T^p\Q) s ® |0) B + Vp\°)s® 
|1>s®|0) b -)■ Pl 1 )^® |0) B + ^|l) s ® |2) B (1) 

Under such a uniquely quantum mechanical noise pro- 
cess, the system S 1 remains in the initial state with the 
probability p of scattering the environment to its excited 
state \1)e if S is in state |0)s and to \2)e if S is in 
state |l)s- h physically describes, for instance, the case 
when a photon scatters randomly in a fibre. The coher- 
ence of S degrades exponentially, and p = 1 — cxp(— Tt), 
where T is the decay rate. When it extends to the case 
of bipartite systems coupling with this environment, the 
decoherence-free states [32| and the phenomenon of en- 
tanglement collapse and revival [33] have been demon- 
strated. 

When an initial state pab = a|$+)($ + | + 
c|vP+)(vP+| +d|*-)(vf-|, with |$±) = 1/V2(|00) ± |H)) 
and 1^) = 1/^2(101) ± ]10>) representing the four Bell 
states and a + b + c + d — 1, evolves in the phase-damping 
channel given by equation ([1]), its off-diagonal elements 
decay exponentially according to the previous analysis. 
Because S(pa) = S(ps) — I, the total correlation is 
calculated as 



T(p 



AB 



2 + ±X 

3=1 



(2) 



where Xj is the four eigenvalues of the final density ma- 
trix. To calculate the classical correlation [(| , the state in 
mode B is projected into \l) — cos0|O) +sin6*e^|l) to get 
the minimal conditional entropy of subsystem A, where 
< 6 < tt and < <f> < 2ir. In the biqubit case, the 
projective measurement is the positive-operator-valued 
measure to maximize C(pab) [34j . which is expressed as 

m 

C{pab) = S(p A ) - min[S(p l A )] 

t),<p 

= ^log 2 (l- ?7 ) + ^log 2 (l+ ? 7), (3) 

S(pa) represents the conditional entropy of A with pro- 
jecting B in \l) and 77 = max{|a|, |/3|, I7I} with a = 
(1 - p)(a - b + c - d), /3 = (1 - p)(c - d - a + b) and 
7 = a + & — c — d. The quantum correlation is therefore 
given by 



Q(pab) = 2 + A 3 lo 82 Aj - C(pab)- (4) 

3=1 



Experimental demonstration of the dynamics of 
correlations. Photon qubits with polarization encoded 
as the information carriers have been widely used to im- 
plement different quantum information processing pro- 
cedures 1]. The coupling between photon polarization 
and frequency modes in a birefringent environment leads 
to dephase by a trace-over frequency [35[, which simu- 
lates the decoherence effect in equation ([T]). Here, we 
encode the horizontal (H) and vertical (V) polarizations 
of a photon as |0)s and \l)g of a qubit and pass the 
polarization-entangled photons through one-sided con- 
trollable quartz plates to investigate the dynamics of dif- 
ferent kinds of correlations. 
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FIG. 1: (Color online). Experimental setup. Ultraviolet (UV) 
pulses are divided into two parts by the ultraviolet polariza- 
tion beam splitter (UV PBS). The relative power between 
them can be changed by a half-wave plate (HWP1). En- 
tangled photon pairs are prepared by spontaneous paramet- 
ric down conversion by pumping the two adjacent crystals 
(/3-barium borate, BBO). Quartz plates (CP) are used to 
compensate the birefringence of down conversion photons, 
both of which then pass through half-wave plates (HWP2 
and HWP4) with the optic axes set at 22.5°. The photon 
in mode A further passes through a half-wave plate (HWP3) 
with the optic axis set as horizontal to prepare the required 
mixed state. The dephasing environment is simulated by 
quartz plates (Q) with thickness L. An unbalanced Mach- 
Zehnder device in the dotted pane with the reflective part 
further passing through two half-wave plates (HWP) is in- 
serted into mode B to prepare another initial mixed state. 
Quarter- wave plates (QWP), half- wave plates (HWP), and 
polarization beam splitters (PBS) are used in each mode to 
set the detecting bases for implementing state tomography. 
Single photon detectors equipped with interference filters (IF) 
are used to detected these two photons. 

Figure 1 shows our experimental setup. Ultravio- 
let pulses are frequency doubled from a mode-locked 
Tksapphire laser with wavelength centred at 780 nm, 
with a 130 fs pulse width and a 76 MHz repetition rate. 
They are distributed into two paths by an ultraviolet 
polarization beam splitter, which transmits 45° linearly 
polarized photons {l/\/2(H + V)) and reflects —45° lin- 
early polarized photons (l/y/2(H — V)). The relative 
power between these two paths can be changed easily by 
a half- wave plate (HWP1), and the time difference be- 
tween them is about 6 ns. They then combine again by a 
beam splitter, and both pump two identically cut type-I 
/3-barium borate crystals, with their optic axes aligned in 
mutually perpendicular planes, to produce polarization- 
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entangled photon pairs [36j| . After compensating the 
birefringence effect between H and V with quartz plates 
(CP), the short pump beam produces the maximally en- 
tangled state |$ + ), whereas the long pump beam pro- 
duces |$~). 

Two HWPs (HWP2 and HWP4), with the optic axes 
set to 22.5°, can change H into l/y/2(H + V) and V into 
1/V2(H - V). The HWP (HWP3) with the optic axis 
set as horizontal in mode A is used to introduce 7r-phase 
between H and V. As our detection traces over the time 
difference between these two production processes, as has 
been shown earlier [37i HI] , the prepared state becomes 



p AB = b\$-)($-\+d\y-)(*-\ ) 



(5) 



where b is determined by the relative power between these 
two pump beams and b+d = 1. Then the photon in mode 
A with frequency distribution f(u>) passes through the 
dephasing channel, which is stimulated by quartz plates 
(Q) with the optic axis set to be horizontal and with 
a thickness of L. Therefore, the decoherence parame- 
ter K = J f(u>) exp(zTo;)daj imposes on the off-diagonal 
elements of equation ([5]), where r = LAn/c, with c rep- 
resenting the vacuum velocity of the photon and An is 
the difference between the indices of refraction of H and 
V. In this case, we get that p = 1 — \k\ in equation ([!]). 

An unbalanced Mach-Zehnder device in the dotted 
pane, which further separates the photon into long and 
short paths, is inserted into mode B to prepare another 
input state. The long path in the dotted pane contains 
a HWP with the optic axis set at 45° and another HWP 
with the optic axis set as horizontal. The time difference 
between these two paths is much larger than the coherent 
time of the photon and less than the coincidence window 
of the logic circuit (at about 3 ns). As a result, the pre- 
pared state becomes 



p = dR\$ + ){$ + \ 
+ 6i?|tf + )(vP + | 



■b(l- R)\<S>-)(<$>- 
d(l -R) | 



(6) 



where R is the effective total reflectivity of the two par- 
tial reflecting mirrors in the dotted pane. The photon 
in mode A then further passes through the dephasing 
channel. Finally, the evolved state is reconstructed by 
tomography. Quarter-wave plates, HWPs and polariza- 
tion beam splitters are used in each mode to set the usual 
16 measurement bases 39]. These two photons are de- 
tected by single-photon detectors equipped with 3 nm 
(full width at half maximum) interference filters. 

Figure 2 shows the correlation dynamics of input state 
§5§ with b = 0.75 in the phase-damping channel. The- 
oretically, the minimum of conditional entropy S(p l A ) in 
mode A is obtained by projecting the photon in mode 
B onto \l) = cos9\H) + sm9e l ^\V), with optimization 
over the angles 9 and </>. In our experiment, we measured 
S(p A ) as a function of 9 with (f> = in different thickness 
of quartz plates, as shown in Fig. 2a. We found that 
the minimum value of S(p A ) is obtained with 9 — 45° 
at about L < I38A and 9 = 0° at about L > 138A 
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FIG. 2: (Color online). The correlation dynamics of input 
state (0 with b = 0.75. (a) The values of conditional entropy 
S(pa), with photons evolving in different thickness of quartz 
plates (L) , as a function of the measurement angle 9 in mode 
B. (b) The dynamics of correlations. Green upward-pointing 
triangles, black squares, red dots, blue stars and magenta 
downward-pointing triangles represent experimental results of 
I, C, Q, En and Rn with the green solid line, black solid line, 
red dashed line, blue dotted line and magenta dotted line 
representing the corresponding theoretical predictions. Non- 
entanglement quantum correlation (T>) is further compared 
with C in the inset (the x axis represents the total thickness 
of quartz plates and the y axis denotes different kinds of cor- 
relations). Purple dots represent the experimental results of 
T> and the purple dotted line is the corresponding theoreti- 
cal prediction. Error bars correspond to counting statistics. 
An = 0.78 pm. 



(An = 0.78 pm is the central wavelength of the photon), 
which agree well with theoretical predictions (solid lines) 
[23}. This implies that there is a sudden change in be- 
havior in the decay rate of C according to equation . 
The dynamics of total correlation (I), classical correla- 
tion (C) and quantum correlation (Q) are shown in Fig. 
2b. We find that C (black squares) decay monotonically 
at L < 138An, and then remains constant at L > 138An, 
which agrees with Fig. 2a. Such immunity of classical 
correlation against decoherence implies an operational 
way to compute classical and quantum correlations [23] . 
In contrast, Q (red dots) behaves in the opposite manner; 
it remains constant at L < 138An and decays monotoni- 
cally at L > 138Ao. C and Q overlap at the specific thick- 
ness of about L — 138 An. It is interesting to observe the 
"decoherence-free" area of quantum correlation, which is 
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due to the equal decay rate of I and C. This phenomenon 
is later shown explicitly in theory by Mazzola et al. [4(| ■ 
It may have important applications for those quantum in- 
formation protocols using only quantum correlation. We 
can see that I (green upward-pointing triangles) decays 
exponentially all the time, which is consistent with the 
continuous dynamics of quantum mutual information ac- 
cording to equation @. The green solid line, black solid 
line and red dashed line are the corresponding theoretical 
predictions of the total, classical and quantum correla- 
tions. Error bars correspond to counting statistics. 

The dynamics of entanglement is also shown in Fig. 
2b, which is both characterized by the entanglement of 
formation ( En) (4lj and the relative entropy of entangle- 
ment (Rn) [42| (the analytical expressions of both char- 
acterizations in our case are given in Methods). The blue 
stars in Fig. 2b are the experimental results of En and 
the blue dotted line is the corresponding theoretical pre- 
diction. The value of entanglement is set to at a thick- 
ness of about L — 173Ao, which shows the phenomenon 
of entanglement sudden death |43j , whereas the magenta 
downward-pointing triangles are the experimental results 
of Rn. Because the final state in our experiment can be 
transformed into a Bell-diagonal form, with a local uni- 
tary operation that does not change its entanglement, 
the analytical expressions of Rn (see Methods) are used 
to give the theoretical prediction, which is represented by 
the magenta dotted line. Although the En is larger than 
the relative Rn at the beginning of evolution, they suffer 
from sudden death at the same thickness, thereby con- 
firming their self-consistence. We can see that quantum 
correlations can be smaller or larger than the En during 
evolution. Specifically, quantum correlation decays ex- 
ponentially, whereas entanglement disappears completely 
at the thickness of L > 173Arj. This confirms the early 
prediction that quantum discord is more robust against 
decoherence than entanglement [ljl [2(| • Because the rel- 
ative Rn is on the equal footing of other measures of 
correlations in the form of entropy [29| . it is not always 
larger than quantum correlations, as verified by our ex- 
perimental results. The inset in Fig. 2b further com- 
pares non- entanglement quantum correlation defined as 
V = Q — Rn [201 with classical correlation. Purple dots 
are the experimental results of T>, with the purple dot- 
ted line representing theoretical prediction. Because of 
the sudden change in behaviour of Q, the decay rate of T> 
also suffers from sudden changes and we find that T> < C. 

Another kind of correlation dynamics is shown in Fig. 
3, in which the sudden change behaviour disappears. The 
input state is obtained by equation ([5]) with b — 0.5 (the 
separated state, Fig. 3, top panel) and b = 1 (the maxi- 
mally entangled state, Fig. 3, bottom panel). The min- 
imum of S(p l A ) is gotten with 9 — 45° in Fig. 3 (top 
panel) and 9 = 0° in Fig. 3 (bottom panel) (</> = 0) [23f . 
The Green upward-pointing triangles, black squares and 
red dots represent the experimental results of I, C and 
Q, respectively. The quantum correlation in Fig. 3 (top 
panel) remains at and the total correlation is equal 



P 



o 
O 




50 100 150 200 250 300 
Total thickness of quartz plates (X ) 



FIG. 3: (Color online). The correlation dynamics of input 
state ((5| with b — 0.5 (top panel) and 6=1 (bottom panel). 
Green upward-pointing triangles, black squares and red dots 
represent the experimental results of I, C and Q, respectively. 
The green solid line, black solid line and red dashed line are 
the corresponding theoretical predictions. En (blue stars) is 
shown for comparison, with the blue dotted line representing 
the corresponding theoretical prediction. Orange diamonds 
represent the experimental results of A and the orange dotted 
line is the theoretical prediction. A < in a and En is set to 
(blue stars are not shown). Rn (magenta downward-pointing 
triangles) is also shown in b with the magenta dotted line 
representing the corresponding theoretical prediction, which 
completely overlaps with quantum correlation. Error bars 
correspond to counting statistics. 



to the classical correlation, both of which decay mono- 
tonically. Figure 3 (bottom panel) shows another case 
in which classical correlation remains at 1 and quantum 
correlation decays exponentially. Compared with the dy- 
namics of entanglement En (blue stars), we find that en- 
tanglement also remains at 0, as shown in Fig. 3 (top 
panel) (experimental results not shown), which is equal 
to Q there, whereas it is larger than Q during evolution, 
as shown in Fig. 3 (bottom panel). Orange diamonds 
represent the value of A (if A > 0, it represents the quan- 
tity of concurrence [44| ; see Methods) , which also decays 
exponentially. The value of En is set to when A < 
according to its definition. The relative Rn (magenta 
downward-pointing triangles) is also shown in Fig. 3. It 
remains constant at 0, as shown in Fig. 3 (top panel) (ex- 
perimental results not shown) and then decays exponen- 
tially, as shown in Fig. 3 (bottom panel). Rn completely 
overlaps with quantum correlation in both cases, and the 
non-entanglement quantum correlation reads at 0. The 
experimental results agree well with the corresponding 
theoretical predictions. 

We further show an interesting case in which Q can 
be larger than C. In this case, the dotted pane is in- 
serted into mode B, as shown in Fig. [TJ and the initial 
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FIG. 4: (Color online). The correlation dynamics of input 
state (|6]) with b — 0.9 and R = 0.9. Green upward-pointing 
triangles, black squares, red dots, blue stars and magenta 
downward-pointing triangles represent the experimental val- 
ues of I, C, Q, En and Rn, respectively. The green solid line, 
black solid line, red dashed line, blue dotted line and magenta 
dotted line are the corresponding theoretical predictions. The 
inset (the x axis represents the total thickness of quartz plates 
and the y axis denotes different kinds of correlations) shows 
the comparison of T> (purple dots) with C. The purple dotted 
line is the theoretical prediction of T>. Error bars correspond 
to counting statistics. 



input state is given by equation We set b = 0.9 and 
R = 0.9 to maximize the value of Q — C. Fig. 4 shows 
our experimental results. The total correlation (green 
upward-pointing triangles) decays exponentially, whereas 
the classical (black squares) and quantum (red dots) cor- 
relations exhibit sudden changes in behaviour in their 
decay rates at about L = 78 Ao. After this special thick- 
ness, the classical correlation remains constant and the 
quantum correlation decays exponentially with a higher 
decay rate. During the evolution, Q is observed to be 
larger than C at the thickness interval from about 50Ao to 
90Ao, within error bars. Such an observation contradicts 
the early conjecture that C > Q [3(| and it is consistent 
with previous theoretical predictions [H, |25|, H^]- The 
dynamics of entanglement, which are characterized both 
by En (blue stars) and Rn (magenta down triangles), are 
also shown in Fig. 4. Both characterizations suffer from 
sudden death at the same thickness of about L — 202Ao 
and the entanglement is set to in the subsequent evo- 
lution. The value of the non-entanglement quantum cor- 
relations V (purple dots), which also exhibit a sudden 
change in the decay rate, is compared with the classical 
correlation in the inset. We find that entanglement and 
non-entanglement quantum correlations are always less 
than the classical correlation in our experiment. 



III. DISCUSSION 

The classical and quantum correlations considered in 
this experiment are defined by the one-partition mea- 
surement of a bipartite system @, 0|- In such cases, an- 
alytical solutions for classical and quantum correlations 
have been obtained for some kinds of states with a high 
symmetry [H, [13, IH1 ■ Although classical and quantum 
correlations may be generally asymmetric according to 
the choice of partition to be measured, for states with 
maximally mixed marginals, they are symmetric under 
the interchange of the measured side [HI . There is also a 
quantifier of classical correlation with measurement over 
both partitions of a bipartite system, which is defined 
as the maximum classical mutual information [t| [Toj . 
For states with maximally mixed marginals, the defini- 
tions of classical and quantum correlations are numeri- 
cally verified to be equal to that defined with one-side 
measurement [2~ij . A thermodynamic approach is also 
developed to quantify correlations [1, [Tl|. In particu- 
lar, the quantum-information deficit is defined to char- 
acterize the difference between the information that can 
be localized with Closed Local Operations and Classical 
Communications and the total information of the state, 
and it quantifies the quantumness of correlations 
We find that the quantum-information deficit is equal to 
quantum discord for Bell-diagonal states. On the basis 
of the idea that classical states can be measured with- 



out disturbance, Luo [13[ use the measurement-induced 
disturbance to characterize classical and quantum cor- 
relations. Because of the symmetric property of Bell- 
diagonal states, these definitions of correlations used in 
our experiment also coincide with that defined by Luo's 
method fl4| . Recently, a method to quantify quantum 
and classical correlations on an equal footing with rela- 
tive entropy has been proposed [ljjj and it is applicable 
to multi-particle systems of arbitrary dimensions. In our 
case, considering the dynamics of Bell-diagonal states, 
the one-side definitions of classical and quantum correla- 
tions [a, [7J are consistent with those defined by the equal 
footing method [l4|. As a result, for the Bell-diagonal 
states used in our experiment, all of the above separa- 
tion methods for classical and quantum correlations are 
coincident. 

Through measuring classical and quantum correlations 
in multipartite systems, Bennett et al. [45| recently in- 
troduced three postulates that should be satisfied by any 
measure or indicator of genuine multipartite correlations. 
They found that the concept of covariance proposed by 
Kaszlikowski et al. 46] to identify the existence of gen- 
uine multipartite correlations does not satisfy two postu- 
lates, and they concluded that it cannot be an indicator 
of the genuine multipartite classical correlations [45j . Be- 
cause the equal footing method mentioned above can be 
directly applied to multipartite systems with a unified 
view of quantum and classical correlations it may 
inspire insight into current debates. 

The relationship between the En and quantum cor- 
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relation is also interesting. In our experiment, the En 
is observed to be even larger than quantum correlation 
in certain areas of evolution. For bipartite systems, it 
has also been shown that quantum mutual information 
considered as the measurement of total correlation can be 
smaller than the En [U |48| , which may occur in the case 
of particles with dimensions larger than five (48j . Such a 
contradiction is due to the different frameworks of their 
definitions. In our case, we can choose the relative Rn 
to measure entanglement, so as to make it on an equal 
footing with other measures of correlations in the form 
of entropy [29| . 

In conclusion, we have shown the sudden change in 
behaviour in the decay rates of different kinds of corre- 
lations in the two-photon system. Our results show that 
classical and quantum correlations can remain unaffected 
under certain decoherence areas, which will play impor- 
tant roles in distinguishing correlations [23| and designing 
robust quantum protocols based only on quantum corre- 
lation. The early conjecture that C > Q [30( is disproved 
by our experimental results. On the other hand, classical 
correlation is always found to be larger than the entangle- 
ment and non-entanglement quantum correlation in the 
experiment. 

IV. METHODS 

Analytical expressions of the En and the rela- 
tive Rn. For biqubit systems, the En can be given by 



the analytical formula f 1 il 



1 -L /l "y~2 

En(p)=H( +V 2 ~ ) (7) 

where H(x) = — x\og 2 x — (1 — cc)log 2 (l — x). T is 
the concurrence, given by T = max{0,A}, where A = 
VxT - y/X2 - y/X3 - y/Xi and Xj are the eigenvalues 
in decreasing order of the matrix p{p2 ® G-i)p*{p2 ® 02) 
with (72 denoting the second Pauli matrix and p* cor- 
responding to the conjugate of p in the canonical basis 
{\HH), \HV), \VH), \VV)}. In addition, the relative Rn 
is defined as the minimal distance between p and a dis- 
entangled state [42}. For the Bell-diagonal state, if the 
four eigenvalues of the density matrix are A* =1 € [0, 1/2], 
the relative Rn is given by [42| 

Rn{p) = 0, (8) 

whereas for any > 1/2 we obtain 

Rn{p) = Xj log 2 Xj + (1 - Xj) log 2 (l - Xj) + 1. (9) 
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